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Abstract 

Networks of weakly nonlinear oscillators are considered with diffusive and time- 
delayed coupling. Averaging theory is used to determine parameter ranges for which 
the network experiences amplitude death, whereby oscillations are quenched and the 
equilibrium solution has a large domain of attraction. The amplitude death is shown 
to be a common phenomenon, which can be observed regardless of the precise nature 
of the nonlinearities and under very general coupling conditions. In addition, when the 
network consists of dissimilar oscillators, there exist parameter values for which only 
parts of the network are suppressed. Sufficient conditions are derived for total and 
partial amplitude death in arbitrary network topologies with general nonlinearities, 
coupling coefficients, and connection delays. 

PACS codes: 02.30.Ks; 84.35. +i; 87.10. +e; 05.45. +b 

Keywords: coupled oscillators; time delay; stability; nem'al networks 

1 Introduction 

The interaction of oscillatory systems arise naturally in many areas of science and has 
intrigued investigators for a long time. The earlier studies go back to the 17th century when 
the Dutch scientist Christiaan Huygens observed that two clocks mounted on a common 
basis tend to synchronize regardless of the initial conditions [1, 2]. In the last several 
decades, models based on coupled oscillators have found diverse applications in physics 
and biology. Among many examples are arrays of Josephson junctions [3,4], semiconductor 
lasers [5,6], relativistic magnetrons [7], chemical reactions [8-11], circadian pacemakers 
[12, 13], intestinal electrical rhythms [14], and a variety of biological processes [15-17]. 
The investigation of these systems has provided striking examples of different types of 
dynamical behavior that can be induced by the presence of coupling. Synchronization is 
probably the most widely studied among these [18-20]. Another interesting dynamics is the 
so-called amplitude (or oscillator) death, whereby individual oscillators cease to oscillate 
when coupled and go to an equilibrium solution instead. One of the first accounts of this 
phenomenon was given in the 19th century by Lord Rayleigh, who noted that when two 
organ pipes stand side by side they "... may almost reduce one another to silence" [21, 
Vol. 2, Chapter XVI]. Later, amplitude death was observed in electronic systems [22,23] 
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and chemical oscillators [24], causing renewed interest in this area. A theoretical analysis 
was given in [25], which considered a pair of linearly coupled oscillators each of which is 
near a Hopf bifurcation. In particular, it was established that a fairly general class of 
such systems can exhibit amplitude death in case the individual oscillator frequencies are 
sufficiently different and the coupling is diffusive. In a similar vein, it was shown that a 
network of limit cycle oscillators can experience amplitude death if there is a distribution 
in frequencies [26-28]. 

A series of more recent studies have demonstrated that time delays in the coupling can 
induce amplitude death even in the case of identical frequencies [29,30]. In these works, a 
pair of identical Stuart-Landau equations are considered under scalar and diffusive coupling. 
A linear stability analysis and numerical simulations are used to show that the presence of 
delay can stabilize the equilibrium solution in this case. A similar conclusion is drawn for an 
all-to-all coupled network of identical Stuart-Landau systems, where the coupling strengths 
and time delays between pairs are also identical. Prom a different perspective, other studies 
have shown the role of feedback delays in suppressing or modifying oscillations in a single 
oscillator [31-33]. Time delays are thus significant in the oscillations and stability of many 
systems, which justifies the efforts devoted to their analysis. 

The purpose of the present paper is to obtain rigorous conditions for amplitude death 
when weakly nonlinear planar systems (such as those near a Hopf bifurcation) are diffusively 
coupled and the coupling involves time delays. We consider arbitrary nonlinearities and 
network structure to prove the generality of amplitude death for the interconnection of such 
systems. Furthermore, we show the existence of parameter ranges for which only parts of 
the network are quenched while the others continue to oscillate, even though they remain 
connected. This phenomenon, which may be termed partial amplitude death, is shown to 
be possible in case when the individual units are sufficiently different. 

This study extends the previous results in several directions. Firstly, we treat a general 
class of systems and derive conditions for delay-induced amplitude death irrespective of the 
particular nonlinearities that may be present. In the absence of delays, general conclusions 
can be drawn by considering normal forms that represent the dynamics on a center manifold 
when each oscillator is near a Hopf bifurcation [25]. However, the analysis is more difficult 
when delays are introduced, and as a result, the study of amplitude death in the delayed 
case so far has been confined to some special equations. For instance, the system studied in 
[29,30] is obtained by adding delays to the center manifold equations considered in [25], and 
as such it does not represent general center manifold dynamics for delay-coupled oscillators. 
Hence, an interesting qiicstion is what conclusions can be drawn in the presence of delays 
when nonlinearities other than the Stuart-Landau system are considered. 

Secondly, we consider very general coupling conditions. We do not restrict to scalar 
coupling and present results for arbitrary diffusion matrices between oscillators. Both the 
diffusion matrices and the delays are allowed to vary from pair to pair within the network. 
In addition, there are no assumptions on either the size or the topology of the network, and 
the results apply to any network layout, including all-to-all, nearest-neighbor, and sparse 
coupling types. 

Finally, instead of a local stability analysis we obtain global results using averaging 
theory. Among other things, this allows the results to be applicable to more general systems 
such as bistable oscillators, where the origin is stable but there are also limit cycle solutions 
which may be annihilated by the effect of coupling. In fact, in such cases the local stability 
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analysis of the origin cannot yield parameter ranges for amplitude death since the origin 
always remains stable. On the other hand, the averaging approach enables one to study 
the (averaged) trajectories inside a large sphere and derive conditions that ensure that 
they all decay to zero. Furthermore, this non-local view is essential for investigating the 
phenomenon partial amplitude death, since here also, the stability type of the origin for 
the overall system remains unchanged as parts of the network collapse to approach the zero 
solution. 

The paper is organized as follows. In the next section a pair of coupled oscillators is 
considered. After a brief discussion of averaging theory for delay-differential equations, the 
averaged equations are obtained, which indicate that qualitatively there are two cases to 
consider. Section |21 treats the first case where the frequencies of individual oscillators are 
close to each other. Here, sufficient conditions are derived for the occurrence of amplitude 
death, and the possibility of death for identical oscillators for almost all positive values of 
the delay is proved. The case of large frequency differences is addressed in Section IH where 
it is shown that it is possible for only one of the oscillators to be suppressed while the 
other one continues to oscillate. Section El extends the results to an arbitrary size network 
with varying coupling coefficients and delays between the oscillators. A general condition is 
derived under which the network experiences amplitude death as a whole or only in parts. 
Results are illustrated on numerical examples. 

For notation, ||-|| is used for the Euclidean norm. If IT is a matrix, then denotes its 
transpose, tr its trace, and Kg = \{K'^ + K) its symmetric part. An identity matrix of 
appropriate size is denoted by /, while J denotes the skew-symmetric matrix 



Furthermore, the derivative of a function / is written as Df. 



2 A pair of coupled oscillators 

We consider oscillators that can be modelled by equations of the form 

Xi{t)=AiXi{t)+£fi{Xi{t)), (1) 

where x, G R^, the matrix Ai G R-^^^ has a pair of purely imaginary eigenvalues iiwi 7^ 0, 
and e is a nonnegative parameter. The function fi includes any nonlinear terms and is 
assumed to satisfy the following hypotheses. 

(H) The function : R^ ^ R^ is and /,(0) = 0. 

It is further assumed that Ai is put into the real Jordan form 

(2) 

with LOi > 0, possibly after a linear coordinate transformation. Equation typically 
arises from a center manifold reduction of higher order systems near a Hopf bifurcation. 



A; = LOi J 





-UJi 



Ui 
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Alternatively, it may be viewed as a perturbation of the harmonic oscillator, which includes 
the classical pendulum, van der Pol, and Duffing equations as special cases. 

Denoting ^i{t) = exp{tAi), a transformation of variables given by Xj = ^i{t)ui puts (CQ) 
into the form 

Ui = e<^-\t)M'^,{t)ui). (3) 

This time-dependent equation can be made autonomous by averaging theory. The averaged 
equation corresponding to Q is defined as 

Ui = efi{ui) , (4) 

where ^ 

/.(n) = lim i / ^-\t)n{Mt)u)dt. (5) 

T^oo 1 Jq 

By the averaging theorem [34], nonzero and hyperbolic equilibria of Q correspond to 
periodic solutions of for all sufficiently small e. If one considers limit cycle oscillators, 
((IJ) will usually have one or more asymptotically stable equilibrium points and the origin will 
be unstable; however, these assumptions are not essential for the purposes of the present 
paper. 

Now suppose a pair of systems of the form (P) are coupled diffusively, 

xiit) = Aixiit) + eh{xi{t)) + £K{x2{t - t) - xi{t)) , . 

X2{t) = A2Xi{t) + e/2(x2(i)) + eK{xi{t - r) - X2{t)) ^ ' 

Here K G R^^^ is a matrix of coupling coefficients and r is a nonnegative number repre- 
senting the transmission delay. The matrix K is scaled by the parameter e to indicate that 
the coupling strength is assumed to be of the same order of magnitude as the attraction 
to the limit cycles of individual oscillators. Following the change of variables Xi = ^i{t)ui, 
((HI) becomes 



iiiit) = e$r'(t)/i($i(t)ni(t)) + e^^\t)K ($2(t - r)n2(t - r) - 

Mt) = e^2 ' W/2(^2(i)^2(t)) + e^^\t)K - r)ni(t - r) - Mt)u2(t)) 



For small £, ui and U2 are slowly changing variables and the method of averaging is appli- 
cable. 

Averaging theory for functional differential equations such as (j7| is in principle similar 
to the theory for ordinary differential equations. To introduce some notation, let y denote 
the pair (^1,^2) G R", with n = 4 in the present case. Then ((T)) is a relation between 
the derivative y(t) and the present and past values of y{t). A solution y{t) of describes 
a trajectory on the state space C = C([— r, 0], R"), the space of continuous functions 
mapping the interval [— t, 0] to R". The trajectory consists of points yt which are defined 
by yt{d) = y{t + 0), 6 £ [— r, 0]. Hence each point yt £ C on the trajectory corresponds to 
a window of the solution y{t) over an interval of length r. The system (j7| is a particular 
case of the more general functional differential equation 

y{t) = £g{t,yt) , (8) 
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where y € R" and g : (RxC) R". For small values e, ^ can be viewed as a perturbation 
of a trivial ordinary differential equation. In fact, if the following limit exists 

1 

g{ip) = hm - / g{t, if) dt , 
then the averaged equation corresponding to (jHI is defined to be 

in which yt is interpreted as a constant element of C. In this way, averaging reduces the 
infinite-dimensional system to an ordinary differential equation on R" [35]. In the 
following, we shall exploit this fact for the analysis of 0. The main averaging result is 
contained in next lemma. 



Lemma 1 Let $j = exp(tAj), where Ai is given by and let K G R . Then 



T 



K^j{-T) if LOi = LOj 

if uJi^ LOj 



where 



K=i(trK)./-itr(JK). J. 



Proof. By (jSJ one has 

= exp(t^i 



cos ijjit sin ijjit 
— sincjjt cosa;,t 



(9) 

(10) 

(11) 



Clearly, ^~^{t) = ^f{t) and $i(t-r) = <^i{t)<^i{-T). Thus, if L denotes the left hand side 

of a, 



L 



lim — 



^J{t)K^j{t)dt 



cl>,(-r). 



(12) 



When oji = loj = uj, the integrand above consists of second order homogenous polynomials 
in the variables sinwt and coswt. Let [kmn] denote the elements of K. Then, using the 
elementary facts that 



1 



T-.00 T 



lim — / sin^ cot dt = lim — / cos^ ujtdt = — 



1 



and 



one can calculate 



lim — / smojt cos ijjt dt = 



lim — 

T->oo T 



^l{t)K^j{t)dt 



ku + k22 ki2 - k2i 
k2i - ki2 ku + k22 
1 



1 

2 

i(tri^)./-ltr(JK). J. 
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So, L = K^j{—T) for the case of equal frequencies. On the other hand if Ui ^ loj, then the 
entries of the matrix ^f[t)K^j{t) are hnear combinations of terms of the form p{ujit)q{ujjt), 
where q can be the sine or the cosine functions. It follows by the theory of almost periodic 
functions [36] that the average in the above sense of such terms is zero whenever tOi ^ loj. 
m 

As the lemma implies, it is necessary to distinguish between two qualitatively different 
cases, determined by the difference \u;i — iV2\ between the frequencies. These cases shall be 
addressed in the next two sections. 



3 Similar oscillators and amplitude death 

We first consider the coupling of two oscillators with similar native frequencies. More 
precisely, suppose that — ui2\ = 0{e). This includes the particular case of identical 
frequencies. In fact, writing Ai = A2+ (^i^), and subsuming the 0{e) difference into the 
definition of the function fi or /2 in @, it is clear that there is no loss of generality in 
taking uji = UJ2, which will be assumed for the rest of this section. 

With uji = UJ2 = oj, applying Lemma ^ to the coupled system © gives the averaged 
equations 

ill = e{fi{ui) - Kui + K<^{-t)u2) ,-.o^ 
U2 = e{f2{u2)-Ku2 + K^{-T)ui) ^ ' 

where $ = $1 = $2 and the fi are defined by 0. We shall show the existence of parameter 
values such that inside an arbitrary ball centered at the origin all solutions {ui{t),U2{t)) of 
(|13j) approach zero as t — > 00. This will imply the amplitude death in the original equations 
^ for sufficiently small e. 

Theorem 1 Suppose K G R^^^ and R > 0, and let the numbers Qi, i = 1,2, be defined by 

= supj^^^j^^ : ||n|| < i?, w / o| . (14) 

// 

tr K — \ (tr K) cos lot — tT:{JK) sin ujt\ > 2 max{gi , ^2} (15) 

then all solutions of satisfy limj_>oo('Ui(t), U2(i)) = whenever ||(ui(0),ti2(0))|| < R. 
On the other hand, if 

tvK <^tT{Dfi{0)+Df2m (16) 

then the zero solution of HL^) is unstable. 

Proof. Under the conditions (H), it is easy to see that the supremum in p4j) is finite. 
In fact, since the integrand in © is periodic in t, 

fi{n) =^J^ ^-\t)f,{^,{t)u)dt (17) 

where T' is the period. Therefore the fi are diff'erentiable and vanish at the origin, and one 
has the finite Taylor series expansion 

Mu) = Df,iO)u + n^iu) (18) 
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where the remainder term TZi is continuous and hm^^o 1 1 ^« (^) 1 1 / 1 1 ^ 1 1 = 0- Hence n^/j (u) / 1 1 n 1 1 ^ 
is bounded inside any ball with the origin removed, so the qi defined by (|14() are finite. Now 
taking the inner product of the first equation in l\YA^ by ui and the second one by U2, and 
adding, one obtains 



^^(11^1 11^ + ||'W2||^) = e {Uifliui) + U2f2{u2) - uJKui - U2KU2) 



+ £ {u{K^{-t)u2 + U2K^{-t)ui) 



(19) 



By (d ujfi{ui) < qi 



Hence, ii v = {ui,U2) denotes the vector formed by concate- 



nating ui and U2, then (|19j) yields 



dt 



1^ < 2ev'^Pv 



(20) 



P 



(21) 



where 

Qil-K K<^>{-t) 
K^{-t) q2l-K 

Clearly, Pv = v'^PgV, where Ps denotes the symmetric part of P. We claim that Pg is 
negative definite if the condition (jlSj) is satisfied. Indeed, it is immediate from ()1U() that 
Kg = ^{tr K)I, and a calculation using (fTT|) shows that 



{K^{-t))s = -[(triT) coswT - tr(Ji^) sinwr]/. 



Hence, Ps has the form 

Ps = 

with 
and 



qil-Ks {K^i-r))^ 



rji = 2qi - tr K, 



1 


r]il 


al 


~ 2 


al 


ml _ 



1,2 



(22) 



(23) 



a = (tr K) cos ujt — tr{JK) sin ujt. 

If (|15j) holds, then rji and r]2 are both negative and r]ir]2 > cx^- Thus, by the well-known 
tests for definiteness, the matrix Pg in 1)231) is negative definite, and the claim is proved. It 
now follows under these conditions that 



v'^Pv 



v^PsV < -A \\vf for all w G R"^ 



where — A < is the largest eigenvalue of Ps- Substituting into ^ 

\\v{t)f < ||t;(0)f exp(-2eAt) for t > 



gives 



which establishes the first statement of the theorem. 

To prove the second statement, consider the linearization of (|13() about zero: 

Df_i{0)-K K^{-t)_ 
K^{-t) Df2{0)-K 



(24) 
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By ()17|) and Lemma ^ 



^r\t)Dfi{0)^i{t)dt 
-^tr(J[D/,;(0)])- J. 



= -tv[DMO)]-I 



(25) 



Since J has zero trace, 



Also, tr = tr K by (|Tn|) 



trL>/i(0) =trZ;/i(0). 
Thus the matrix in (|24|) has trace equal to 



(26) 



tr(D/i(0)+Z)/2(0)-2i^) 



which is positive if (|16() holds, in which case at least one of its eigenvalues has a positive 
real part. ■ 

Remark. The quantities qi defined in (|14j) satisfy the lower bounds 



which provides a connection between the stability condition (jlSj) and the instability 
condition ()16() . The inequality H27|) follows by taking limits in (jlHj) and using H25|) to 
obtain 



relates qi to the local stability of the origin for the individual uncoupled oscillator. 
For instance, if the linear part of the oscillator has the eigenvalues a it iuj, then 
gj > Q by (|27|) . In particular, > in case the oscillations arise from a supercritical 
Hopf bifurcation. On the other hand, if the origin is asymptotically stable then 
trD/j(0) < 0, and qi may be negative in a sufficiently small neighborhood of the 
origin, but not necessarily in larger neighborhoods (see Example 1 below). To avoid 
trivial cases from the point of study of the death of limit cycles, it will be assumed 
in this paper that the qi are positive in the region of interest. 

In essence, Theorem ^ states that in case of equal frequencies, amplitude death occurs 
whenever the trace of the matrix K is sufficiently large, provided that sinwr ^ 0. By the 
averaging theory, for a set of parameter values satisfying (P3|) there exists eo > such that 
the original coupled system © undergoes amplitude death if e E (0, eo)- Condition ((T3|) 
is graphically depicted on the parameter plane in Figure H It should be noted that the 
graph is not uniform in e. For instance, it may be necessary to decrease e as ti K increases 
in order to observe amplitude death. This should be expected since the foregoing analysis 
assumes that the coupling coefficients eK are 0(e) terms. Nevertheless, the value of Eq can 
be chosen to be constant over any compact subset of the parameter space. The extent of 
the death region in the figure suggests that amplitude death is quite a robust dynamical 
behavior for diffusively coupled systems. 



qi>-tv[DfiiO)] 



(27) 
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Figure 1: Regions marked "D" are parameter values for which amphtude death occurs, and 
Q = maxjgi, 52}- The sohd Unes correspond to the case when the matrix K is symmetric. 
For a nonsymmetric K, the death regions are warped to the side, as shown by dotted hues. 



Some general qualitative characteristics of the region of amplitude death can be deduced 
from Figure E For instance, the death region is a disconnected set, formed of "islands of 
death" as mentioned in [30]. Furthermore, it is bounded away from the horizontal axis, 
implying that there is a threshold of coupling strength below which amplitude death is 
not possible. Strictly speaking, these properties do not follow from H15() . which is only a 
sufficient condition for amplitude death. Nevertheless, they can be established directly by 
a linear stability analysis of the origin. Thus, the statement that tr should be sufficiently 
large for amplitude death follows from the instability condition (|16() . Similarly, it can be 
shown that amplitude does not occur if ojt = nvr for some integer n. Indeed, in this case 
$(— r) = (—1)"'/ by Hill) , and the linear variational equation (|24() takes the form 



Z)/(0) - K 
{-lYK 



Dm - K 



where for simplicity it is assumed that the oscillators are identical, and / 
Letting C denote the matrix on the right side, we note that 



/l = /2. 



If A is an eigenvalue of Df{0) corresponding to an eigenvector p, then 



I 




I 


_ (-1)"/ _ 


p = 





p 
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i.e., A is also an eigenvalue of C. From (|25|) it is seen that itr[D/(0)]±iitr(J[Z)/(0)]) are 
the eigenvalues of Df{0), and thus of C. Consequently, the origin is unstable if tic[Df(0)] > 
0, e.g. when the oscillations result from a supercritical Hopf bifurcation, as noted in the 
Remark above. This proves that amplitude death is not possible for identical oscillators if 
UT = mr, n £ Z; in particular death does not occur for r = 0, which agrees with earlier 
results obtained for the undelayed case [25]. It should be pointed out that a disconnected 
death region in the parameter space arises under the assumptions that e is small and the 
delays are discrete. For large values of e single or multiple death islands (or none) may 
be possible, as was shown for the delayed Stuart-Landau equations in [30]. On the other 
hand, recent results indicate that regions of amplitude death typically grow and merge if 
the delays are distributed over an interval rather than concentrated at a single point [37]. 

The global nature of Theorem^makes it applicable also to cases where a linear stability 
analysis may be inconclusive — for instance, to bistable oscillators, as the next example 
illustrates. 

Example 1 Consider the following variants of van der Pol oscillators which are coupled 
through their velocities: 



zi + eii(l - (1 + p'^)zl + y^zf) + zi = e/3(i2(t - vr/2) - ii(t)) 
+ ez2{l - (1 + p^)zl + Ip^zi) + za = ep{h{t - vr/2) - iajt)) 



For /9 = and e < the left hand sides of H28() describe the usual van der Pol 
oscillators. The scalar (3 denotes the coupling strength. Putting (|28j) into the vector 
form © by letting Xj = (zj, ij), it is seen that 



K 





(3 



thus tr K = p. In the absence of coupling, the averaged equations have the form 
with 



f^(u) = --i^l-^^^\\ur + !^\\urj-u, i = l,2. (29) 

The roots u of fi satisfy ||m|| = 0,2, or 2/ p. For definiteness we take p = 2. Then for 
P = and £ small and positive, each (uncoupled) oscillator has two attractors: the 
zero solution and a periodic solution with amplitude near 2. In addition, there is an 
unstable periodic solution whose amplitude is near 1. Now by (|29|) . 

1 ' - ^ " '|2 I l II ||4 \ II ||2 



Mu) = -- j^l - - llnf + J llnf J • (30) 
from which, by simple calculus on the quantity in parenthesis, it follows that 

u^Mu)<-^\\u\\'^, ueR^, i = l,2 

so qi = q2 = 9/32. Theorem ^ then implies that for /? > 9/16 and sufficiently small 
positive £ the coupled system has no periodic solutions and the origin is the only 
attractor. The agreement with the numerical solutions of H28() is shown in Figured 
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Figure 2: Death of the periodic solutions of coupled bistable oscillators as the coupling 
strength (3 is varied. The zero solution remains asymptotically stable throughout the pa- 
rameter scale. 



for several values of e. Here the vertical scale is the amplitude (in the variables zi 
and Z2) of the attracting periodic solution, which exists for the uncoupled system, 
and thus also for small values of (3. As the coupling strength (3 increases past a 
critical value, the periodic solutions are annihilated. The critical coupling strength 
approaches 9/16 « 0.56 as e is made smaller, as predicted by averaging theory. A 
closer inspection of the averaged equations reveals the mechanism through which 
amplitude death occurs. With r = 7r/2, Lemma Ogives K = ^(31 and K^{—t) = 
— ^j3J. The latter being an antisymmetric matrix, the last two terms in (|19|) vanish, 
leaving 

^|(hif + WM') = e {ulh{u^) + ulUm) - ^/3(||nif + \\u2f) 

Since the oscillators are identical and the coupling is symmetric, one could look for 
equal amplitude solutions, say ||ui||^ = ||ii2||^ = r'^ ■ Using (|3()j) in the above equation 
then gives 

r = -\{fi+l-\r' + \r')r'MLh{r). (31) 

The equilibrium solutions of H31() are r = and r = | it ^\/9 ~ 16^, whenever the 
radicands are nonnegative. Thus for — 1 < /? < 9/16 there are two positive equilibria 
corresponding to the unstable and stable periodic solutions of the oscillators. As 
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-2 

0.5 1 1.5 2 

r 

Figure 3: The shape of the right hand side of for various values of [i. The roots 
correspond to the amphtudes of the stable (•) and unstable (o) periodic solutions. As 
(3 is increased the two positive roots coalesce and disappear. The critical value for this 
saddle-node bifurcation is /? = 9/16. 

13 — > 9/16^, these two equilibria approach each other and disappear in a saddle-node 
bifurcation. Correspondingly, for (3 > 9/16 there are no periodic solutions. The situa- 
tion is graphically shown in FigureEl It is seen that the sufficient condition /3 > 9/16 
given by Theorem ^ is also necessary for amplitude death in the present case. Fur- 
thermore, throughout the whole range of (3 values the origin remains asymptotically 
stable. Hence, the local stability analysis of the origin reveals no information about 
the death of the periodic solutions in this example. 

4 Dissimilar oscillators and partial amplitude death 

Suppose now that \uji — W2I = C(l)) i-e-, the frequencies differ by a nonzero constant which 
does not depend on the parameter e. Applying Lemma ^ to ((T)), the averaged equations 
are found to be 

ill = e{fi{ui) - Kui) 
U2 = e{f2{u2) - KU2) 

Note that the equations are decoupled and do not depend on the delay r. In analogy with 
Theorem n the following result can be stated. 

Theorem 2 Let R > 0, and qi he defined by il4\ )- If 

tTK>2qi (33) 
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then all solutions of the i-th equation in satisfy \\m.t^ooUi{t) = whenever ||ni(0)|| < 
R. On the other hand, if 

ivK <tTDfi{Q) (34) 
then the zero solution of i-th equation in \3^) is unstable. 

Proof. Taking the inner product of Ui and the i-th equation in ((SH), and using ((TU]l 
and (|14|) . one obtains 

< [qil - K)ui 
= euf{qil - K)sUi 
= ^{Qi - ^trK) \\ui\\'^ 

from which the first statement of the theorem follows. On the other hand, the linearized 
equation 

ill = e{Dfi{0) - K)ui 

is unstable if tr(L'/j(0) — K) is positive. But the latter quantity is equal to tr(L'/j(0) — K) 
as observed in the proof of Theorem which establishes the second statement. ■ 

Similar to Theorem ^ the above result implies that the pair of oscillators undergo 
amplitude death if tr K is sufficiently large, but this time regardless of the value of the delay. 
Comparison of the conditions (|15|) and (|33|) shows that amplitude death becomes more 
likely when the coupled oscillators have a frequency difference, which agrees with previous 
studies of amplitude death caused by frequency distributions [27, 28] . More importantly, 
in this case the conditions for amplitude death are separate for each equation. Thus, it is 
possible that one of the oscillators is suppressed independently of the other one although 
they remain connected, leading to a partial amplitude death in the coupled system. The 
following example illustrates this behavior. 

Example 2 Consider a pair of coupled van der Pol oscillators 

zi + e{zl - ai)zi + ujIzi = e[5{z2{t - 1) - zi{t)) 
Z2 + e{zl- a2)z2 + ojlz2 = e(3{zi{t-l) - Z2{t)) 

where > and /3 G R. As in the previous example, letting Xi = {uJiZi,Zi) shows 
that P is equal to the trace of the matrix in ©. In the averaged equations H32|) 
one has 

fiiu) = ^(^ai-^\\u\\'^^ -u, i = l,2 (36) 

which confirms that in the absence of coupling the i-th oscillator has an attracting 
limit cycle solution Zi = cos Uit + 0(e). Furthermore, yields 

T I f \ ^ W l|2 
U fi[u) < y \\U\\ . 

Hence, with qi = ai/2 aii application of Theorem [21 gives the condition (3 > Ui for the 
death of the i-th oscillator in the coupled system. In particular, if 02 < /3 < ai then 
only the second oscillator will be suppressed. The results of the numerical solution 
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Figure 4: Behavior of coupled van der Pol oscillators whose limit cycles have different 
amplitudes and frequencies, as the coupling strength /3 is varied. Partial amplitude death 
occurs for 0.5 < /3 < 1, where only one of the oscillators is quenched while the other one 
continues to exhibit high-amplitude oscillations. The region /3 > 1 corresponds to the usual 
amplitude death of the whole system. Note that the vertical scale is logarithmic. 

of (|35|) is summarized in Figure 0J where the parameter values are ai = 1, a2 = 0.5, 
uji = 0.1, UJ2 = 1, and e = 0.1. It is seen that amplitude death occurs for /? > 1, 
where both oscillators are quenched. On the other hand, a partial death occurs for 
0.5 < P < 1, where only the second oscillator is damped while the first one continues 
to exhibit high-amplitude oscillations. For this parameter range the amplitude of the 
second oscillator is almost (but not exactly) zero; in fact, the amplitude approaches 
zero as e — > O"*". Indeed, the averaged equations (|32|) in this case are 



which have an attracting fixed point of the form (ui, 0) with ||ni|| = lyj a\ — f3 when 
02 < (3 < ai. Thus for small e the amplitudes of zi and Z2 in the original system 
(jSSI) are 2-y/l — P + 0(e) and 0{e), respectively. It should be remarked that, for this 
range of P, although the zero solution is asymptotically stable in the equation for U2 
it is not stable for the full system of equations (|37j) . As the origin remains unstable 
throughout the range < /? < 1, a linear stability analysis alone is insufficient to 
reveal the occurrence of the partial amplitude death of the coupled system. Also note 
that the actual value of the delay is not significant in this behavior, and one obtains 
pictures similar to Figure |3 even for zero delay. 




(37) 
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It is rather interesting to observe partial death under diffusive couphng, since one 
intuitively expects this type of connection to act towards reducing the differences between 
the oscillators, when, for instance, X is a positive scalar. For the undelayed case, results 
related to this behavior was given in [38], which essentially follow by the observation that 
weakly connected systems near a Hopf bifurcation can be decoupled in the 0{e) terms if 
the frequencies are different, a fact which has been applied to the study of neural networks 
[39,40]. Here we see that the same conclusion also holds in the presence of arbitrary 
coupling delays. However, it should be noted that higher-order terms in e may still be 
coupled and dependent on r. This can be significant at large values of e; for instance, the 
results of [30] indicate that in such a case amplitude death may depend on the value of 
r, especially when the frequency difference between the oscillators is small. In closing we 
mention that a form of partial death was also observed numerically in undelayed networks 
of relaxation oscillators [41]. Since Theorem |21 holds under quite general conditions, partial 
amplitude death can be expected to be among the common behaviors of many practical 
networks formed by the interaction of dissimilar units. 



5 Network of oscillators 

The results of the previous sections will now be extended to a general network of N coupled 
oscillators modelled in the form 

N 

Xiit) = AiXi{t) + efi{xi{t)) + e ^ Kij{xj{t - nj) - Xi{t)) , i = 1, 2, . . . , iV. (38) 

Here Xi G R^, Ai = tOiJ, fi satisfies the hypothesis (H), and e > 0. The nonnegative quan- 
tity Tij is the transmission delay and Kij G R^^^ are the coupling or diffusion coefficients 
between the i-th and the j-th oscillator. To prevent self-coupling Ka is set to zero for all i. 
It is assumed that the coupling between pairs is symmetric, that is, Kij = Kji and Tij = Tji 
for all On the other hand, both the coupling coefficients and the delays are allowed 
to vary from pair to pair. This is an especially realistic assumption if the delays arise as a 
result of spatial separation of the oscillators. 

With the change of variables Xi = ^i{t)ui = exp{tAi)ui, (|38|) is transformed to 

N N 

iiiit) =e<^r\t)fii^i{t)ui{t))-e^^r\t)Kij<^iit)ui + e^^^\t)^^^ 

j=i i=i 

(39) 

for i = 1, 2, . . . , A^. If r = maxjrjj \ i, j = 1, . . . , N}, then (|39j) has the form (jH)) where the 
right-hand side a function on RxC([— r, 0], R^^), and can be averaged as before. To this 
end, suppose there are M distinct frequencies loi, . . . ,ljm among the oscillators, and let In 
denote the set of indices of those oscillators having frequency w„, for n = 1, . . . , M. More 
precisely, j € In if and only if the matrix Aj has eigenvalues HcOn- It follows by Lemma H 
that the average of the terms of the form ^^^{t)Kij^j{t — Tij) in H39() will be zero whenever 
i £ In and j ^ In- Application of the lemma gives the averaged equations 

/ ^ _ \ 
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for n = 1, . . . , M, where 



Kij = ^(tr Kij) • / - ^ tr{JKij) ■ J. (41) 



The averaged equations consist of M decoupled sets of equations. The results of the 
previous sections suggest that each set can undergo amplitude death independently of the 
others. This observation is made precise by the following theorem. 

Theorem 3 Suppose K £ R^^^ and R > 0, and let qi, i G X„, be defined by \14^ . Let Vn 
denote the vector formed by the concatenation of Ui for i £ Zn- If 

N 

'^tr Kij - ^ \{ti Kij) cosLUnTij - tv{JKij) sinuJnTijl > 2qi for every i G J„ (42) 

then all solutions of the n-th set of equations i40j) satisfy liuit^oo Vnit) = whenever 
||t'n(0)|| < R- On the other hand, if 

i^tv Kij - tr D f,{0)] <0 (43) 
*ex„ \j=i J 

then the zero solution of i4(J[ ) is unstable. 

Proof. Taking the inner product of the i-th equation in H4U|1 with Ui and adding over 
i £ In gives 

Let i? > 0. By ((EI) and (dU, 



dt II "I - ^ \ II 'II 2 , 



provided maxjgj^ ||nj|| < R. Let S denote the last summation in i.e., 

5 = ^ ^ ujK,j<^j{-Tij)Uj 

= Y.Y. \nJ[K^,^j{-r^,) + ^]{-n^)KJ^Uj 
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where the last hne follows since Kij = Kji and Tij = Tji for all and = for i,j G X„ 
by our assumptions. Hence from (|22() . 

where 

Oij = (tr Kij) cos LOnTij — tl{JKij) sinLJnTij. (46) 

From the identity 

< \\ui ± UjW'^ = + \\uj\\'^ ± 2ufuj 

follows the estimate 

II l|2 I II l|2 
\Ui Uj\ < 

Using in 1)45^ and noting that aij = aji gives 
Substitution into yields 



^ij I W'^il 



dt 



< eA ^ ||nj| 



2 



where we have defined 

N 



A = max 2qi — tr i^j^ + | 
i=i i6i„ 



Thus 

lbn(t)f < |b„(0)f exp(eAt). 

If condition ()42() holds, then A is negative, and the first statement of the theorem is proved. 
To prove the second statement, consider the linearization of (|4U() about zero: 

/ N _ \ 

= £ Dfi{0)ui - ^KijUi + ^ Kij^j{-Tij)uj , i G Tn. 
V i=i iex„ / 



Recalling that Ka = 0, and using (|26|) and (|4T|) . the divergence of the vector field defined 
by the right hand side is calculated as 



N \ / N 
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which is positive if H43|) holds. In this case, any sufficiently small neighborhood of the origin 
expands in volume under the evolution of the differential equations H4U|) . which proves the 
instability of the zero solution. ■ 

In view of the above theorem it is worthwhile to examine the death of a given network 
as additional units are connected to it. Suppose that the network consists of A'' oscillators 
with identical (or 0(e) close) frequencies, and consider the condition (|42|) for amplitude 
death. When an additional oscillator with a different frequency is connected, the first sum 
in H42|) changes by tr(ifj^Ar+i) while the other terms remain the same. Thus if tr(Ki^isf^i) 
is sufficiently large for i = 1, . . . .A^, then the original set of A'^ oscillators will experience 
amplitude death after the new connection. Conversely, if tr(i^j^jv4.i) is negative for some i, 
then it is possible for the new connection to "awaken" the network from a previous death 
state, by satisfying the instability condition (jlH)) . 

When the newly added oscillator has the same frequency as the rest, the coupling 
delay becomes important. For simplicity let us assume that all oscillators and the coupling 
conditions between them are identical, so the indices can be dropped. The death condition 
(|42|) then has the form 

trK- |(trK)cosa;T-tr(Ji^)sinLjr| > 2q/N. (47) 

Assuming q > 0, this condition cannot be satisfied for any A^ if either sinwr = or 
tr K < 0. However, if sinwr ^ and K is symmetric with positive trace, the network 
will undergo amplitude death for all sufficiently large network sizes A^. For nonsymmetric 
matrices the same conclusion holds whenever the left hand side of (|47)) is positive, i.e., if 
tr K is sufficiently larger than \tr{JK)\. However, for highly nonsymmetric matrices the 
presence of the connection delay may make amplitude death impossible, as the left hand 
side of (|^7|l becomes negative. 

Although physically meaningful in most models, the assumption that there are no self- 
connections in the network (|38|) is not an essential requirement for Theorem El In fact, an 
investigation of the first part of the proof shows that H42|) remains a sufficient condition for 
amplitude death even when the matrices Ka are not zero. This observation points to an 
interesting connection between the present paper and previous works on feedback control 
of oscillations. For if one considers the system 

x{t) - Ax{t) - ef{x{t)) = eK{x{t - r) - x{t)) (48) 

as a special case of the network ()38() consisting of a single oscillator. Theorem El gives the 
sufficient condition 

tr K — \ (tr K) cos ujt — tr{JK) sin ujt\ > 2q 

for its stability. Since the latter inequality can be satisfied if tr K is large and sinwr ^ 0, one 
concludes that for almost all positive values of the delay r and sufficiently small e > 0, the 
oscillator given by the left hand side of can be stabilized by a linear delayed feedback 
of the form K{x{t — r) — x{t)), regardless of a precise knowledge about the nonlinearity 
/. This is in agreement with related results which were obtained for similar systems in the 
context of delayed-feedback control of oscillatory behavior [31-33,42]. 

The possibility that the network may be divided into clusters which separately experi- 
ence amplitude death entails the formation of spatial patterns consisting of oscillating and 
suppressed units. Such pattern formation is best studied under more specific assumptions 
on the network structure; therefore, it is not treated further in this paper. 
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6 Conclusion 



We have presented a mathematical analysis of the phenomenon of amplitude death for a 
general class of weakly nonlinear systems which are diffusively coupled with time delays. 
Sufficient conditions are derived for the coupled network to exhibit amplitude death. It 
is shown that when the individual oscillators are sufficiently different, it is possible for 
only parts of the network to be suppressed, while the rest continues to oscillate. The size 
of the network may enhance or inhibit amplitude death depending on the nature of the 
coupling matrices, where both the magnitude and the sign of the trace of the matrix plays 
an important role. The quantitative conditions presented are rigorously justified for weakly 
nonlinear oscillators. In particular, they arc significant when considering networks where 
each individual unit is near a Hopf bifurcation. On the other hand, numerical simulations 
suggest that in many cases the qualitative aspects of the results remain valid for stronger 
nonlinearities as well. In view of the few assumptions required on the network structure 
and the individual units, the generality of the results imply that the death of the network, 
either in part or as a whole, may be among the common and robust dynamical behaviors 
in many practical systems. This should be significant in studies of interacting periodic 
processes, such as physical, biological, and financial cycles. 
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